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Transversity quark distributions in a covariant quark-diquark model 
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Transversity quark light-cone momentum distributions are calculated for the nucleon. We utilize 
a modified Nambu-Jona-Lasinio model in which confinement is simulated by eliminating unphysical 
thresholds for nucleon decay into quarks. The nucleon bound state is obtained by solving the 
relativistic Faddeev equation in the quark-diquark approximation, where both scalar and axial- 
vector diquark channels are included. Particular attention is paid to comparing our results with the 
recent experimental extraction of the transversity distributions by Anselmino et al. We also compare 
our transversity results with earlier spin-independent and helicity quark distributions calculated in 
the same approach. 
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I. INTRODUCTION 

For perfectly collinear quarks, at leading twist, the nu- 
cleon has three independent quark distribution functions 
for each quark flavour: the unpolarized distributions, 
q{x), the helicity distributions, /S.q(x), and the transver- 
sity distributions, Axqix). Knowledge of each distribu- 
tion is of equal importance if one is to have a robust de- 
scription of nucleon structure. The unpolarized and he- 
licity distributions have been extensively studied, both 
experimentally and theoretically, for many years [l|, 
However interest in the transversity distributions is rela- 
tively recent 

The transversity distributions are associated with 
quark-nucleon forward Compton amplitudes, where both 
the quark and nucleon hclicitics arc flipped. Hence these 
distributions are chiral-odd. The electroweak and strong 
interactions conserve chirality, so the transversity distri- 
butions must couple to another chiral-odd quantity in 
scattering cross-sections. This is not possible in inclu- 
sive deep inelastic scattering (DIS) However, for cer- 
tain semi- inclusive DIS [1, IH, Q and Drell-Yan processes 
B B 0, IF] the transversity distributions do appear at 
leading twist in the cross-section. For semi-inclusive DIS 
the transverse quark distributions couple to particular 
chiral-odd fragmentation functions (Collins functions), 
whereas in Drell-Yan dilepton production they appear 
with either a transverse quark or anti-quark distribution 
from the partner hadron. 

The transversity distributions are particularly interest- 
ing for several reasons, for example: 

• the moments of the transversity valence distribu- 
tions are related to the nucleon tensor charge @; 

• the quark bilinears associated with transversity are 
odd under charge conjugation and hence the quark- 
anti-quark sea does not contribute Q; 

• helicity conservation at the quark-gluon vertex 
prevents mixing between the quark and gluon 



transversity distributions under QCD evolution 

• gluon transversity distributions, Axgix), only ex- 
ist for targets with J ^ 1, because measurement of 
a gluon transversity distribution requires that the 
target change helicity by two units of angular mo- 
mentum and this is not possible for spin-i targets 

These results imply that the transversity distributions 
are valence quark dominated and evolve as non-singlets 
under DGLAP evolution, where the angular momen- 
tum generated by the DGLAP kernels is not shared 
between the quark and gluon sectors. These features 
of the transversity distributions make them particularly 
amenable to a quark model treatment. 

In this paper we calculate the transversity distribu- 
tions for the proton using a Nambu-Jona-Lasinio model 
[H, [ll] . This model is attractive because it is covariant 
and has a transparent description of spontaneous chiral 
symmetry breaking. Confinement - in the sense that 
there exists no threshold for nucleon decay into quarks 
- is also implemented via the regularization procedure 
[13, [iBl . We construct the nucleon as a bound state so- 
lution of the relativistic Faddeev eq uation [T6| . in the 
quark-diquark approximation p^l . Il8l |. where both scalar 
and axial-vector diquark channels are included. We com- 
pare our transversity results to spin-independent and he- 
licity quark distributions calculated in the same approach 
p^ . Particular attention is also paid to a comparison of 
our transversity results with the recent, and to date the 
only, experimental extraction of the transversity distri- 
butions by Anselmino et al., presented in Rcf. [20|. 



II. TRANSVERSITY QUARK DISTRIBUTIONS 

The transversity distributions render a probability in- 
terpretation analogous to the other two leading twist 
distributions. In a transversely polarized hadron they 
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represent the number density of quarks in an eigenstate 
of the transverse Pauh-Lubanski operator, i^_l75, with 
eigenvalue minus the number density of quarks with 
eigenvalue — i , fsj that is 

Arqix) ^ q^{x) ~ qi{x). (1) 

In an helicity basis the helicity distributions are expressed 



Aq{x) = q+{x) - q-{x), 



(2) 



where q+{x) is the number density of quarks with he- 
licity parallel to the hadron helicity and q~-{x) is the 
quark number density with helicity anti-aligned. The 
spin-independent distributions in each basis are given by 



q{x) = (Ji (x) + qi (x) ^ q+ {x) + g_ {x) . 



(3) 



The leading twist quark distribution functions are de- 
fined by light-cone Fourier transforms of connected ma- 
trix elements of particular quark field bilinears. For ex- 
ample, the twist-2 transversity distribution is defined by 



^Tq{x)=p+ / |1 



(p,s|^^(0)7+7'75V',(r)b,s)c, (4) 

where tpq is a quark field of flavour q and x is the Bjorken 
scaling variable.^ In Eq. ([4]) the target polarization is 
in the cc-direction, with the z-direction defined by the 
photon 3-momentum. 

The evaluation of the quark distributions is facilitated 
by expressing Eq. (|4]) in the form [2l|, 

Arqix) ^-ij -0^6 (^x - Tr [7+7^75 M(p, k)] , 

. . (5) 

where M{p, k) is the quark two-point function in a nu- 
cleon. The quark distributions can then be expressed 
in terms of Feynman diagrams for any model where the 
nucleon is represented by a bound state of quarks. The 
diagrams we consider are given in Fig. [T] In our pure va- 
lence quark model there should also be a third diagram, 
the so-called quark exchange term , however this dia- 
gram does not contribute within the static approximation 
used here 

In the Feynman diagrams of Fig. [1] the single line rep- 
resents a constituent quark propagator and the double 
line a diquark t-matrix. The diagram on the left is re- 
ferred to as the quark diagram and on the right we have 



^ The formal expressions for the spin-independent and heUcity dis- 
tributions can be obtained from Eq. Q via the operator replace- 
ments 7~''7^75 — > 7"*" and 7~''7^75 — > 7~''75, respectively. The 
nucleon state is normalized according to the non-covariant light- 
cone normalization, namely {p,s\ipy^^'^ipu -I- V'd7^'0d|P> = 3. 
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FIG. 1: Feynman diagrams representing the transversity 
quark distributions in the nucleon, needed in the evaluation 
of Eq. ([5}. The single line represents the quark propagator 
and the double line the diquark f-matrix. The shaded oval 
denotes the quark-diquark vertex function and the operator 



insertion has the form 7^7^75 5^ a 



i(l±r. 



the diquark diagram, where we include both scalar and 
axial- vector diquarks. Separating the isospin coefficients, 
the u- and d-quark transversity distributions can be ex- 
pressed as 
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^ATf^^n)/Ni^) + l^TQN{^) 
+ I ATQn)/Ni^) + ^At/,7c)/^(x), (6) 

+ I ^T^DyNi^) - ^Arf^l^y^ix), (7) 

where each term represents a particular Feynman dia- 
gram in Fig.[TJ The superscripts s, a and m refer to the 
scalar, axial-vector or mixing terms, respectively. The 
subscript q/N implies a quark diagram and q{D)/N a 
diquark diagram. Because the scalar diquark has spin 
zero, we have At f^(^jjyj^(x) = and hence the transver- 
sity of the d-quark arises exclusively from the axial- vector 
and the mixing terms. 

Importantly, in this covariant framework, the Ward 
identities corresponding to number and momentum con- 
servation are satisfied, guaranteeing the validity of the 
baryon number and momentum sum rules [23l . from 
the outset. 



III. THE NUCLEON IN THE NJL MODEL 

The NJL model which we employ to determine the 
transversity distributions is discussed in Ref. p^ . where 
it was used to determine the spin-independent and helic- 
ity quark distribution functions. Therefore we give only 
a brief outline of the model in this section and refer the 
reader to Ref. for further details. 

We represent the nucleon as a quark-diquark bound 
state, where both the scalar and axial-vector diquark 
channels are included. The nucleon t-matrix is obtained 
as the solution of the relativistic Faddccv equation 



T ^ Z + ZWnT 



THat Z. 



(8) 
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In the static ap proximation the quark exchange kernel, 
Z, becomes [1^ 



M VV37; 



'5 7^^ 



\fil^ 75 



(9) 



where M is the constituent quark mass. With this 
approximation, Hat of Eq. ([5]) effectively becomes the 
quark-diquark bubble graph, given by 



where 



r(g) 



T.(<Z) 

Tr(<z) 



(10) 



(11) 



The quantities Tgifi) and T^"'((j') are the solutions to the 
NJL scalar and axial- vector diquark Bethe-Salpeter equa- 
tions, respectively [Ulli,!!!], and S{k) in Eq. ^ is the 
usual fermion propagator for a constituent quark. 

In the non-covariant light-cone normalization used al- 
ready in Eq. the quark-diquark vertex function, r^r, 
is defined by the behaviour of T near the pole 



Tat Fat 



(12) 



where = ^ 



is the nucleon's light-cone energy. The 
vertex functions satisfy homogeneous Faddeev equations 
of the form 



N i- N, 



and T M ~ T N Hn Z 



(13) 



As the NJL model is non-renormalizable, a regular- 
ization prescription must be specified to fully define the 
theory. We choose the proper-time regularization scheme 
[11, Ea, HE where loop integrals of products of 

propagators are evaluated by introducing Feynman pa- 
rameters, Wick rotating and making the denominator 
replacement 



Xn 



{n-iy. yi/(Lc/v) 



(14) 



where hm and huv are, respectively, the infrared and 
ultraviolet cutoffs. The former has the effect of eliminat- 
ing unphysical thresholds for hadron decay into quarks, 
hence simulating an important aspect of quark confine- 
ment The proper-time scheme also preserves the 
gauge invariance and Poincare covariance of the theory. 




FIG. 2: Model results for the triplet of twist-2 valence 
up quark distributions, at Qq = 0.16 GeV'^. The spin- 
independent and helicity distributions are taken from earlier 
work presented in Ref. [19t |. 



Ga [iMl • The infrared or confinement scale is expected to 
be of order Lqcd and we choose Ljr = 240 MeV, how- 
ever our results exhibit only a minor dependence on hiR 
when it is varied between 200—300 MeV. The parameters 
TO, Ijuv and the qq coupling in the pion channel, G^r, are 
determined by requiring M = 400 MeV via the gap equa- 
tion, /tt = 93 MeV from the one loop pion decay diagram 
and TOtt = 140 MeV from the pole of the qq i-matrix in the 
pion channel. This gives to = 16.4 MeV, Ljjv = 645 MeV 
and Gtt — 19.04 GeV~^. The qq couplings in the scalar, 
Gs, and axial- vector, Ga, diquark channels are deter- 
mined by reproducing the nucleon mass, = 940 MeV, 
as the solution of Eq. and satisfying the Bjorken 
sum rule within our model, with gA = 1.267. We obtain 
Gs = 7.49 GeV-2 and Ga = 2.80 GeV'^. With these 
model parameters the diquark masses are Ms = 687 MeV 
and Ma = 1027 MeV. 

Inherent in most model determinations of the quark 
distributions is the absence of an explicit scale. Con- 
sequently the model scale, Qq, must be determined via 
comparison with empirical results. We do this by op- 
timizing Qq such that the spin-independent valence up 
quark distribution, w^(x) = u{x) — u{x), best reproduces 
the empirical parametrization after NLO Q^ evolution.^ 
We find Qg = 0.16 Ge V^ |19l. which is typical of valence 
dominated models [1^ l26l. l30l| . 

In Figs. [H and [3] we illustrate our results for the 
transversity valence u- and d-quark distributions, respec- 
tively. Included in these figures are our results for the 
spin-independent and helicity valence distributions taken 
from Ref. ^1^. We find a transversity u-quark distribu- 
tion slightly larger than the corresponding helicity dis- 



IV. RESULTS 



The six parameters of the model are: the regulariza- 
tion parameters, Lm and hjjv, the current quark mass, 
TO, and three quark-quark coupling constants G^, Gs and 



^ The DGLAP evolution of our results for the three twist-2 quark 
distributions is performed using the appropriate computer pro- 

For the DGLAP parameters we 
: 3 and Aqcd = 250 MeV. 



gram from Refs. [3l|, l32l . |32 
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FIG. 3: Model results for the triplet of twist-2 valence 
down quark distributions, at Qq — 0.16 GeV'^. The spin- 
independent and helicity results are taken from earlier work 
presented in Ref. p^ . 



tribution for all x. The transversity and helicity d-quark 
distributions are very similar for x > 0.35, however at 
smaller x the transversity distribution is noticeably sup- 
pressed. The difference between the transversity and 
helicity valence distributions results purely from rela- 
tivistic effects, since in the non-relativistic limit, where 
boosts and rotations commute, there can be no preferen- 
tial polarization direction and therefore the helicity and 
transversity distributions must be equal. This is true 
only at the model scale as the helicity and transversity 
distributions evolve differently under DGLAP evolution. 

The first experimental extraction of the transversity 
distributions was achieved only recently, and is pub- 
lished in Ref. [l^l- The authors combined semi- inclusive 
DIS data from HERMES and COMPASS with e+e" 
annihilation data from BELLE to simultaneously ex- 
tract the transversity distributions and Collins functions. 
Their results for the transversity u- and d-quark distri- 
butions at = 2.4 GeV^ are presented in Fig. [Has the 
shaded regions, which represent a one-sigma confidence 
interval. Included in this figure are our valence results 
at — 2.4 GeV^ and the empirical Soffer bound at 
the same scale obtained from an evolution of the GRV 
parametrizations given in Refs. [13, UHl- A direct com- 
parison between our results and those of Ref. [l^] should 
be valid for x > 0.2 where transversity anti-quark dis- 
tributions are expected to be small. We find that our 
results lie slightly outside the one-sigma error bounds of 
the empirical parametrizations for x > 0.3. On the same 
figure we illustrate our results for the helicity distribu- 
tions at = 2.4 GeV^. At this scale our helicity and 
transversity distributions remain very similar for x > 0.4, 
however the differing evolution has resulted in a sub- 
stantial suppression for the transversity distributions at 
smaller a:, when compared with the results at the model 
scale. The helicity distributions given in Fig. 3] are in 
excellent agreement with the empirical parametrizations 



FIG. 4: The shaded areas are the empirical results of Ref. [201 ] . 
with a one-sigma confidence interval. The dot-dashed line is 
the GRV Soffer bound [s^, HI] and the solid lines are our 
results for the transversity distributions. The dotted curves 
are our helicity distributions taken from Ref. [l9| ]. All results 
are at = 2.4 GeV^ 



The first moments of the transversity valence distribu- 
tions are related to the nucleon's isovector tensor charge, 
9T, via [| 



/ dx [ATUy{x, Q' 
Jo 



)~ATd,{x,Q')]^gT{Q'). (15) 



The nucleon's isoscalar tensor charge, gj^, is defined as 
the sum of the valence transversity moments. For these 
moments we obtain Atu^ = 1.04 and Axd^ = —0.24, 
giving a nucleon isovector tensor charge of gx — 1.28 and 
a isoscalar charge of = 0.80, at — 0.16 GeV^. At 
the GRV scale of = 0.4 GeV^ we obtain Atu^ = 0.69, 
Ard,, = -0.16, gr = 0.85 and g^ = 0.53.3 j|- jg ^gjj 
known that the first moments of the helicity distributions 
give the quark spin content of the nucleon. However, the 
first moments of the transversity distributions are not 
equivalent to the quark spin in the transverse direction, 
since the expectation value of 7^75 is not equal to that of 
Y,± — 7^7^75. That is, the isoscalar tensor charge cannot 
be interpreted as a transverse spin sum Q . 

The paper of Anselmino et al. does not give explicit val- 
ues for the transversity moments, however if we integrate 
their parametrizations and accurately determine the er- 
rors from Fig. 7 in Ref. [13], we find Atu = 0.46^^5:28 and 

= 0.4 GeV^.* Using these mo- 
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n iq+0.30 



at 



^ These values are obtained using the NLO result [I 



1- — 



ATg(Qo)- (16) 



4867r 

* These values are obtained from the FIT-I results given in Ta- 



0.2 0.4 0.6 0.8 1 

X 

FIG. 5: Quark distribution ratios, in each case the dotted line 
is the result before QCD evolution and the solid line is the 
resuh for = 5GeV^ 
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FIG. 6: Quark distribution ratios, in each case the dotted line 
is the result before QCD evolution and the solid line is the 
resuh for = 5GeV^ 



ments their tensor charges are gr ~ 0.65 and gj^ ~ 0.27, 
assuming Atu — Ard ~ 0. 

The similarity between the transversity and helicity 
distributions, which we illustrate in Fig. 31 is a com- 
mon theme for model calculations - see Refs. IH, 113, 
m, m, SH. This shouW come as no surprise since 
the constituent quark masses are typically rather large, 
350 — 450 MeV, in these calculations and therefore enor- 
mous relativistic effects would be unexpected. For the 
moment, we have no more than a hint that this may 
disagree with the experimental extraction of Ansclmino 
et ai, who find a u-quark transversity distribution that 
may in fact be far smaller than its helicity counter- 
part. This difference is emphasized if we study the mo- 
ments. The empirical values for the helicity moments are 
Au = 0.851±0.075 and Ad = -0.415±0.124 42], for the 
valence distributions the results are Au^ = 0.926 ±0.071 
and Ad.„ = -0.341±0.123 [13]. Therefore Atu extracted 
from the fit of Ref. [lO] is potentially far smaller than 
the M-quark helicity moment. If this difference survives 
it would indicate a difficulty for relativistic model ap- 
proaches based on the concept of a constituent quark 
mass. 

Within our NJL model calculation we find Auy — 0.97 
and Ady = —0.30,^ which differ only slightly from the 
transversity results given earlier. The nine models dis- 
cussed in the transversity review by Barone et al. Q all 
find Atu > Au and Axd ^ Ad at the model scale. 



ble I of Ref. [2Cll |. The FIT-II results give even smaller moments, 
namely AyM = 0.40 and A^d = —0.16. The errors quoted in the 
text are first obtained at = 2.4 GeV^, then the NLO result 
expressed in Eq. llbt is used to devolve to the parametrization 
scale, = 0.4 GeV^. We find that this procedure works very 
well for the central value, which can be determined at both 
scales. 

^ Note, the helicity moments given here differ slightly from those 
in Ref. [iSlI as in this work the infrared cutoff is slightly different. 



This consensus amongst the various model approaches, 
together with our results, contrasted with the poten- 
tially small empirical w-quark transversity parametriza- 
tion given in Ref. [lO], presents an interesting problem. 
We therefore eagerly await new experiment data so that 
the empirical transversity distributions can be better con- 
strained. 

The Soffer inequality [ii], q{x) + Aq{x) > 2 \ATqix)\, 
and the more familiar positivity conditions shed little 
light on the size of the it-quark transversity distribution, 
simply that |ATw(a;)| ^ u{x). However, the differing 
signs of the spin-averaged and helicity d-quark distribu- 
tions implies that for all values of x where Ad{x) ^ the 
transversity d-quark distribution must satisfy |ATd(a;)| ^ 
^d{x). For the d-quark moments the Soffer inequality im- 
plies lAydl ^ |, provided Ad ^ 0. In fact, if Ad ^ — ^, 
which appears true empirically [sl, H^, then the 
transversity moment must satisfy lA^dl ^ Our re- 
sults, presented in Figs. [5] and [31 satisfy all positivity 
constraints, including the Soffer bound. 

The relative importance of each diagram in Fig. [1] to 
the quark distributions is easily illustrated if we con- 
sider their contributions to the transversity moments. 
For the -u-quark distributions we have At/^/jv = 0.75, 
At/„V = -0.016, At/:(^)/^ - 0, At/:(^)/^ = 0.090 
and At/^^-i/jv ~ 0.22, where the isospin factors in 
Eq. ([6]) have been incorporated into the numerical val- 
ues. Results for the d-quark are obtained by dividing 
out the M-quark isospin factors and multiplying by the 
d-quark factors from Eq. ([7]). The dominant terms for 
Atu{x) are clearly the scalar-quark and diquark-mixing 
diagrams. The strength of d-quark distribution arises al- 
most exclusively from the diquark-mixing term. Similar 
results are also obtained for the helicity distributions |1_3] . 
This clearly illustrates the importance of including axial- 
vector diquark correlations in any quark-diquark picture 
of the nucleon. 

In Figs. [5l and [H] we illustrate various ratios of the 
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transversity distributions with the unpolarized and hehc- 
ity distributions. In each case the dotted line is the result 
at the model scale, = Qqi ^^'^ solid line is the ratio 
evolved to = 5 GeV^. In Fig. Owe present our results 
for the single flavour ratios of transversity to unpolarized 
quark distributions. Both the up and down quark ratios 
are finite in the x 1 limit at the model scale. However, 
as all transversity moments decrease with increasing Q^, 
we expect the ratios to decrease after QCD evolution, 
which is evident for the = 5 GeV^ results.® 

Ratios of transversity to helicity distributions are pre- 
sented in Fig. [SI In the large x limit we find that the u- 
quark ratio approaches 1.1 and the d-quark ratio becomes 
quite large. The d-quark result is surprising. It may 
be an inherent feature of the model or possibly a conse- 
quence of the static approximation. The quark exchange 
diagram, absent in this approximation, contributes dom- 
inantly in the scalar diquark channel, where a d-quark 
must be exchanged. Therefore, including this contribu- 
tion may significantly alter the ratio Axd{x) / Ad{x) at 
large x. This possibility illustrates the importance of 
going beyond the static approximation in future work. 
The final result illustrated in Fig. [S] is the mixed flavour 
transversity ratio, Axdv^x)/ Atu^^x), which is small for 
all X and tends to zero as x approaches one. 

To aid comparison between our calculations and fu- 
ture investigations of the transversity distributions, we 
provide a parametrization of our u- and d-quark valence 
distributions at the GRV scale, = 0.4 GeV^. The 
general form of the parametrization is 

xArqvix) = rjx°'{l - xf {I + + px). (17) 

For the w-quark distribution we find 

T] = 1.5025, a = 1.1586, (3 = 3.9940, 

7 = 113.65, (5 = 11.150, p= 12.214, (18) 

and a fit to the d-quark distribution gives 

77 = -2.7990, a = 1.5952, /? = 5.7531, 

7 = 45.424, 5 = 8.0163, p = 2.2065. (19) 

These numbers produce an extremely good fit to our re- 
sults, however one should not read too much into their 
values as a similar can be achieved with other sets 
of parameters. The transversity parametrization used in 



The non-conservation of the nucleon tensor charge results from 
a non-zero anomalous dimension in Mellin space, Ax-lgqin), at 
leading order. The sign of the transversity anomalous dimen- 
sion at both the leading and next-to-leading order is negative, 
and therefore the tensor charge drops with increasing Q^. More- 
over, since Ar7gq(n) < A7qq(n) for all moments n, the transver- 
sity moments decrease more rapidly than the helicity moments. 
This has important consequences for the QCD evolution of the 
transversity distributions. 



Ref. [lOl has the form 

xArq^ix) = ^N;[x"{l-xf 

^—L^L_[q{x) + Aq{x)], (20) 

where q{x) and Aq(x) are the GRV parametrizations 
given in Refs. (3^. l35j. Using this parametrization, a 
fit to our distributions at = 0.4 GeV^ gives 

iVj = 1.13, iVj = -1.46, a = 1.48, /3 = 0.70, 

(21) 

which also furnishes a good representation of our results. 
These values are very different from those obtained in 
Ref. [1^, refiecting the fact that our distributions are 
considerably larger. 

V. CONCLUSION 

Using a covariant quark-diquark model for the nucleon, 
which includes both the scalar and axial-vector diquark 
channels, we calculated the transversity quark distribu- 
tions for the proton. A key feature of this framework is 
that it produces quark distributions that have the cor- 
rect support and obey the number and momentum sum 
rules. This model also eliminates unphysical thresholds 
for nucleon decay into quarks and hence incorporates im- 
portant aspects of confinement. 

We found that at the model scale the transversity 
and helicity distributions are very similar in magnitude, 
with AtUv^x) being slightly larger than Auv{x) for all 
X, whereas ATdy{x) becomes smaller than Adv{x) for 
X < 0.35. These results satisfy the Soffer inequal- 
ity. For the first moments of the distributions we find 
AxUy — 1.04 and Ardy — —0.24, giving a nucleon isovec- 
tor tensor charge of gr = 1.28 at the model scale. This 
result is very near the empirical value for the axial charge, 
gA — 1.267, which we use as a constraint to help de- 
termine the model parameters. Axial-vector diquarks 
were found to play a pivotal role, essentially through 
the scalar-axial- vector mixing term included here for the 
first time in a transversity calculation. This diagram al- 
most exclusively gives rise to the d-quark distribution 
and provides about 20% of the strength for the u-quark 
transversity result. Future work in this direction will be 
to calculate fragmentation functions, in particular the 
Collins functions, so that a direct comparison between 
semi-inclusive DIS experiments and theory is possible. 

We compared our transversity results with the re- 
cent experimental extraction of Anselmino et al, and in 
each case found that our distributions lie slightly out- 
side the one-sigma confidence interval of their empirical 
parametrizations, in the valence quark region. The po- 
tentially small magnitude of the empirical Atu{x) distri- 
bution, if confirmed by future experimental data, would 
indicate very large relativistic effects for the quarks in the 
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nucleon. This would be very difficult to explain within 
relativistic models utilizing the concept of a constituent 
quark mass. 
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